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Abstract 


In the past two decades, Global Navigation Satellite System-Reflectometry (GNSS-R) has emerged as 
a new remote sensing technique for soil moisture monitoring. Some experiments showed that the 
antenna of V polarization is more favorable to receive the reflected signals, and the interference 
pattern technique (IPT) was used for soil moisture and other geophysical parameters retrieval. 
Meanwhile, the lower satellite elevation angles are most impacted by the multipath. However, 
electromagnetic theoretical properties are not clear for GNSS-R soil moisture retrieval. In this paper, 
the Advanced Integral Equation Model (AIEM) is employed using the wave synthesis technique to 
simulate different polarimetric scattering at the specular directions. Results show when the incidence 
angles are larger than 70°, scattering at RR polarization (The transmitted signal is Right Hand Circular 
Polarization (RHCP), while the received one is also RHCP) is larger than that one at LR polarization 
(The transmitted signal is RHCP, while the received one is Left Hand Circular Polarization (LHCP)), 
while scattering at LR polarization is larger than that at RR polarization for the other incident angles 

(1°~70°) . There is an apparent dip for VV and VR scattering due to the Brewster angle, which will 
result in the notch in the final receiving power and this phenomenon can be used for soil moisture 
retrieval or vegetation corrections. The volumetric soil moisture (vms) effects on their scattering are 
also presented. The larger soil moisture will result in lower scattering at RR polarization, and this is 


very different from the scattering of the other polarizations. It is interesting to note that the surface 


correlation function only affects the amplitudes of the scattering coefficients at much less level, but it 


has no effects on the angular trends of RR and LR polarizations. 
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1. Introduction 


Soil moisture is the link between surface water and groundwater, which plays an important role 
in terrestrial ecosystems and water cycle systems. As one of the determinants in the status of the 
terrestrial hydrosphere, soil moisture largely controls the evapotranspiration, water migration and 
carbon cycle on the land surface. Soil moisture is also important for understanding and forecasting 
climate change, solving global water cycle, water resources management, and basin hydrological 
model. Meanwhile, it is also one of the necessary conditions for monitoring crop growth, droughts and 
floods. 

Soil moisture has strong spatial and temporal variation heterogeneity. The traditional ground and 
meteorological monitoring are difficult to meet the requirements. The development of satellite remote 
sensing technology provides a new method for monitoring. In the past two decades, GNSS-R has 
emerged as a new promising remote sensing technique due to its unique advantages [1]. Its 
applications extend from meso-scale ocean remote sensing to the remote sensing of snow,ice and land 
surface parameters [2-4]. As for the land scenario, the most important application is soil moisture 
monitoring. GNSS-R is an effective complement to existing radars and radiometers, and can be a 
potential efficient technique for the validation of the L- band soil moisture satellites, such as SMOS 
(Soil Moisture Satellite) [5] and SMAP (The Soil Moisture Active Passive (SMAP) mission) [6]. 
However, as an emerging means of Earth observation, there are many open problems that need to be 
solved, such as the polarizations properties of the reflected signals. 

Polarizations are one of the most important properties of the electromagnetic waves. The 
polarization of a plane wave describes the shape and direction of the space trajectory formed by the 
end of the electric field vector as a function of time. However, different from the traditional linear 
polarizations of the microwave remote sensing, the transmitted polarization of GNSS (Global 


Navigation Satellite System) is RHCP (Right Hand Circular Polarization). When the GNSS signals 


reflect from the surface, their polarizations will change. Kavak et al. [7] pointed out that the reflected 
signals consisted two parts, the circular polarization and the perpendicular polarization, which are 
dependent on the elevation angles. Several experiments have paid attention to the polarization study of 
GNSS reflected signals. In paper [7], the modified GPS receiver was used, and a horizontal-positioned 
antenna of LHCP (Left Hand Circular Polarization) polarization was used to collect the coherent 
energy of the direct signals and the reflected signals. In the SMEX 02/03(Soil Moisture Experiment 
2002/2003) experiments [8], the antennas of the Global Position System-Reflectometry (GPS-R) 
receiver was different. An airborne GPS-R receiver, known as DMR (Delay Doppler Maps Receiver), 
was used with two antennas. The standard zenith RHCP antenna was used to collect the direct signal, 
and the nadir LHCP antenna was for the collections of the surface-reflected signals. Meanwhile, the 
incidence angles were from 15°~35°[9]. 

As early as in 2000, Zavorotny and Voronovich have used the GO (Geophysical Optics Model) s- 
cattering model [2] to study the reflections polarizations of soil moisture [10]. They simulated the 
sensivity of polarization ratios(LR/RR and HR/VR) to the soil moisture at larger incident angles (60 ° 
and 70°) and concluded theoretically that the polarization ratios can be used for soil moisture 
measurement, here LR ,RR, HR and VR present the polarization of transmitted signal as RHCP, while 
the received polarizations are LHCP, RHCP, horizontal (H) and vertical (V), respectively. In the BAO- 
Tower (Boulder Atmospheric Observatory-tower) experiment [11], a multi-polarization receiver 
provided by the NASA Langley Research Center was used. In order to collect the surface reflected 
signals, a low-gain LHCP and four high-gain antennas (V, H, LHCP and RHCP polarizations) were 
used. The incidence angle was 35° and the incident azimuthal angle was 245°. However, the BAO- 
Tower experiment did not support their original theoretical hypothesis [10]. 

At present, the interference pattern technique (IPT) was used for soil moisture and vegetation 
retrieval [12]. An antenna of V polarization was employed, since the angular information would be 
masked due to the Brewster’s angle if an antenna of LHCP was used [13]. LEiIMON (Land Monitoring 
with Navigation Signals) was a long-term GNSS-R experimental campaign, which was funded by 


European Space Agency (ESA) [14]. The experiment was conducted in 2009 using a SAM sensor, 


antenna of RHCP polarization was used to collect the direct signals and antenna of RHCP or LHCP 
polarization was used to collect the reflected signals. The elevation angles were from 40° to 80°. 
However, their theoretical and experimental results indicated that both the reflectivity of RR and LR 
was increasing with the increase of the incidence angle (0~80°), which was very different from the 
commonly used Fresnel theory of circular polarization in GNSS-R study [2,4]. Different from the 
specially designed GNSS-R receiver, the multipath information of the geodetic-quality GPS receivers 
could be used to remotely sense the soil moisture in the top 5cm of soil [15~17]. The in-situ 
measurements showed that the elevation angels with 5°~30° were most impacted by the 
multipath[15~17]. 

In general, most of the GNSS-R soil moisture monitoring is based on the experimental study [8, 9, 
11, 12, 14-17]. It is commonly used the simple combination of geometry to proceed from linear 
polarization to circular polarization [2, 4, 14~17]. RR polarization reflection increases with the 
incidence angle increases, while LR reflection decreases. The trend is very different from the 
LEiMON experimental report [14]. It should be highlighted that only the amplitudes changes are 
considered [2, 4, 14~17], while phase differences are not considered as calculations from linear to 
circular [18]. However, with the development of GNSS-R, more attentions should be paid to 
theoretical models since they can assist in the data interpretation, sensitivity analysis and experiments 
design, and others [19]. In this paper, we focus on the polarization (especially circular polarization) 
dependence of bare soil scattering using an AIEM model [20,21], wave synthesis technique is used for 
arbitrary polarizations calculation in Section 2. Reflectivity of XR polarization and linear polarizations 
versus incidence angels are presented in Section 3, XR polarization means that the transmitted signal 
is RHCP polarization and any polarization that is circularly polarized (RHCP and LHCP) or linearly 
polarized (V, vertical or H, horizontal) is received. From the Mueller matrix, the effects of volumetric 
soil moisture and surface roughness for different polarizations are evaluated in Section 4. Finally, a 


conclusion is given in Section 5. 


2. AIEM model 


To overcome the ionospheric effects, the transmitted signals are RHCP, while the polarizations of 


the transmitted signals changed after reflected from the soil surface. As shown in Section 1, since the 
polarizations of the reflected signals are different, the linear (Vertical and Horizontal) and circular 
(RHCP and LHCP) antennas are used to receive the reflected signals. Therefore, different 
combinations of polarizations need further study in theory. In this section, the commonly used random 
rough surface bistatic scattering model of AIEM has been employed. By using the method of wave 
synthesis technique, the model can calculate the bistatic scattering coefficient of arbitrary polarization 
combination (both linear and circular). In this way, the developed bistatic scattering model can be used 
to study the full polarization characteristics of soil moisture. 

The traditional electromagnetic scattering models from random rough surfaces [19] are the 
Kirchhoff model (KA) model and the Small Perturbation Model (SPM) model. The former is suitable 
for surfaces with small curvatures or high frequency limit, while the latter is applicable to slightly 
rough surfaces or low frequency. The popular Integral Equation Model (IEM) is developed to bridge 
the gap between the SPM and KA. However, its several assumptions limit the accuracy, especially for 
bistatic scattering. The subsequently developed AIEM model is an improvement of the IEM model [20, 
22]. Here we focus on the polarization study using the developed surface scattering model for GNSS- 
R applications. It is pointed out that the received power is mainly coming from the first Fresnel zone 


[23], and therefore only the coherent scattering is considered in our scattering properties study. 


The bistatic scattering coefficient is related to the modified Mueller matrix as shown below [19]: 
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(1) 


where g, 


m is the bistatic scattering coefficient, subscripts q and p are the polarizations of the received 


and transmitted signals, respectively, A is the illuminated area and the modified Mueller matrix M 


can be written in terms of the elements of the scattering matrix S, 
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Please read the reference[22] for more information of the AIEM model. 


According to wave synthesis technique [24], the scattering coefficient for a general transmitted 


and received antenna polarization can be written as, 
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(3) 
where subscript q and p are the polarizations of the received and transmitted signals, respectively, and 
I is the normalized Stokes vectors with unit amplitude and Q is a rotating matrix. The normalized 
Stokes vectors can be expressed in terms of the ellipticity angle x and orientation angle y . 
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(4) 
The first element of I describes the intensity of the electromagnetic wave in the horizontal direction 


and the vertical direction, the second one is the intensity difference between the horizontal and vertical 


directions of the wave vector, the third one is the linear polarization intensity at azimuth angle of +45° 
in the polarization ellipse and the fourth one indicates the degree of circular polarization. In this way, 
the phase information in scattering matrix S can be calculated through I. 


For like polarization, , and for a cross polarization, — 
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Any polarizations of transmitted and received signals can be obtained by setting the 


ellipticity X and orientation angles % , accordingly. For GNSS-R, it is necessary to study the 
scattering characteristics of Cie: Subscript XR polarization means that the transmitted signal is 


RHCP polarization and any polarization that is circularly polarized (RHCP and LHCP) or linearly 
polarized (V, vertical or H, horizontal) is received. 


3. Reflectivity of different polarizations versus incidence angels 


When a plane wave is transmitted through a plane boundary between the non-destructive media, 
there is an angle at which the vertical polarization broadcasts a full projection, which becomes the 
Brewster angle. The IPT method expressed in Section 1 uses the position and number information of 


the notches in the direct signal and the reflected signals coherent waveform to perform the geophysical 


parameter retrieval. In fact, the notch in the waveform is the result of Brewster angle. Due to its 
particularity, Brewster's angle is a phenomenon that needs attention as for GNSS-R remote sensing 
inversion. Table 1 presents the corresponding dielectric constants and the Brewster angles for different 
soil moisture. The vms is a way of expressing soil moisture, that is, the percentage of soil moisture 
volume. This kind of representation can reflect the degree of water filling of soil voids, and can 
calculate the three ratios of soil solid, liquid and gas. Fig.1. shows the Fresnel reflectivity of rv and rh 
polarization for different vms, here r means Fresnel reflectivity, while subscript v and h means the 
polarization state, v stands for vertical polarization and h stands for horizontal. The corresponding 
Brewster angles in Table 1 are apparent at V polarization (left panel). As we can see if the vms is 
certain, the rv decreases with the incidence angle increase (before the Brewster angle), and the rv 
increases after the Brewster angle. As for the different vms, the rv is larger if the vms is larger (before 
59°), but the rv is smaller if the vms is larger (after 82°). However, the rn is larger if the vms is larger. 

The real parts of the rv and m (Re[tv],Re[rn]) play an important role in the calculations of XR and 
linear polarization scattering coefficients. So their changes with the incidence angles are presented in 
Fig.2. Before hitting the Brewster angle, the Re[ry] is positive, while it becomes negative if the 
incident angle is lower than the Brewster angle, the Re[rv] increases with the incidence angles increase. 
While Re[rn] is negative for the whole angles with Re[rn] decrease as vms increase. 

Fig.3 is reproduced from the paper [25] and it shows the scattering coefficients of linear 
polarization and XR polarization versus the incidence angels. As we can see, for smaller incidence 


angles (Pe 0 <70 ), the scattering of LR polarization is larger than that of RR polarization 
(RR<LR).Here LR and RR polarization mean the transmitted signals are RHCP, while the received 
ones are LHCP and RHCP, respectively. However, for larger incidence angles (70° < 0 <85), the 


trend is just opposite, that is to say the scattering coefficients of RR polarization are larger than the 
one of LR polarization. The scattering coefficients of RR polarization is larger than that of VR 


polarization for larger incident angles (70 <0<85°)} here VR polarization means the transmitted 


signal is RHCP, while the received one is vertical polarization. The scattering coefficients of RR 


polarization increase as the incidence angles increase (1° < @ < 55 ), but for the other incident angles, 


the scattering coefficient of RR polarization decrease. However, the scattering coefficients of LR 
polarization decrease as the incident angle increases. The trends of scattering coefficient of VV, VR, 
HH and HR versus incident angles are the same. Here VR and HR means the transmitted signals are 
RHCP polarization, while the received signals are vertical and horizontal polarization, respectively. 
VV and HH mean the same polarization of vertical and horizontal, respectively. The reason is that 
they are only related to M (1,1) and M (2,2), respectively (Eq2,5). As the incident angles increase, the 
scattering coefficients of HH and HR polarizations decrease. Scattering coefficients of VV and VR 
polarizations have dips at about 70° because of the Brewster angle defined in Eq.8[19], the scattering 
coefficients decrease at first(before the Brewster angle) and then increase with the incidence angle 


increase(after the Brewster angle). 


4. Effects of volumetric soil moisture and surface roughness on polarizations 


4.1 Effects of volumetric soil moisture on polarizations 


Fig.4 gives the linear and circular polarization scattering coefficients at different soil moisture 
(vms) versus the incident angles. The scattering trends of HR polarization are very similar to that of 
HH polarization, while VR polarization is very similar to VV polarization as shown in Fig.1. The dips 
at VV and VR polarization are related to the vms, and this is due to the Brewster angles (Table 1). If 
the vms is larger, the position of dip is at the position of larger incident angle. It should be noted that 
the scattering coefficients of RR polarization are larger if the vms is smaller. For the other 
polarizations (LR,HR,HH), larger vms corresponds to larger scattering coefficients. But for VV and 
VR polarization, if the vms is larger the scattering coefficients are larger if the incident angles are 
lower than the Brewster angles. While if the incident angles are larger than the Brewster angles, the 
lower vms corresponds to larger scattering coefficients. 

4.2. Effects of surface roughness on different polarizations 
Since the soil moisture affects the GNSS reflection signal, the geometric characteristics of the 


soil, that is, the surface roughness will also affect the GNSS reflection signal. Therefore, in the soil 


moisture inversion, how to remove the surface roughness is a difficult problem. The surface roughness 
can be represented by two parameters, the root mean square height (s) and the surface correlation 
length (1), which respectively define the surface roughness from the vertical and horizontal scales. 
Another way to describe the degree of random surface roughness is the surface autocorrelation 
function, which can be described by a Gaussian autocorrelation function or an exponential 
autocorrelation function. 

In this section, surface roughness effects on the final scattering at RR and LR polarization are 
presented. Fig.5 shows the roughness effects on different correlation functions. As we can see, if root 
square height height(s) is fixed, larger correlation length (1) correspond to larger scattering coefficients 
(both RR and LR); If / is fixed, larger s corresponding to larger scattering coefficients (both RR and 
LR), k is the wave number. 

The scattering of exponential correlation function are larger than the one of Gaussian function 
(both RR and LR) (Fig. 6). It is interesting to note that unlike at backscattering direction, the surface 
correlation function only affects the amplitudes of the scattering coefficients, while the angular trend 


is at much less level. 
5. Conclusions 


As for GNSS-R remote sensing, the polarization of transmitted signals are RHCP, which is 
commonly thought that the received power is coming from the first Fresnel zone. Here the coherent 
scattering of different polarization is simulated by the modified AIEM model using the wave synthesis 
technique. Results show that the scattering at RR polarization is lower than that at LR polarization as 
for incident angles lower than 70° , when the incident angles are between 70° and 80° , the 
scattering at LR polarization is larger than that of RR polarization. As for larger incident angles, there 
are dips at VV and VR polarization due to the Brewster angles, which are related to the soil moisture. 
Therefore, the antenna of v polarization is suggested to be used because the information of dips caused 
by the Brewster angle at different soil moisture are useful for soil moisture monitoring, and the 
scattering dynamic ranges of VR polarization are larger than those of the other polarizations. 
Meanwhile, VR polarization will provide efficient information for vegetation corrections. It should be 


noted that larger soil moisture corresponds to lower scattering at RR polarization, while the trend is 


just the opposite as for the other XR (LR, HR, VR, HH, VV) polarizations. Meanwhile, it is 
interesting to note that unlike at backscattering direction, the surface correlation function only affects 


the amplitudes of the scattering coefficients, but it has no effects on the angular trend. 
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Appendix A 
AIEM: Advanced Integral Equation Model 


DDM: Delay-Doppler Maps 

DMR: Delay Doppler Maps Receiver 

GNSS: Global Navigation Satellite Systys 

GNSS-R: GNSS-Reflectometry 

GO: Geophysical Optics Model 

GPS: Global Positioning System 

HR, VR, +45R and -45R polarizations: The polarization of the transmitted signals are RHCP with 
the received ones H( horizontal), V(vertical) +45 (+45° linear) and -45(-45° linear) 
polarizations, respectively. 

IEM: the Integral Equation Model 

IPT: the Interference Pattern Technique 

KA: the Kirchhoff model 

LEiMON: Land Monitoring with Navigation Signals 

LHCP: Left Hand Circular Polarization 

LR: The polarization of the transmitted signal is RHCP, while the received one is LHCP 
RHCP: Right Hand Circular Polarization 

RR: The polarization of the transmitted signal is RHCP, while the received one is also RHCP; 
SMAP: Soil Moisture the Active and Passive Mission 

SMEX: Soil Moisture Experiments 

SMIGOL: Soil Moisture Interference pattern GNSS Observations at L-band Reflectometer 
SMOS: Soil Moisture and Ocean Salinity Mission 

SPM: the Small Perturbation Model 

Z-V model: Zavorotny-Voronovich model 


Table.1 Dielectric constants and the Brewster angles for different vms 


vms Dielectric constant Brewster angle 
0.01 2.814-0.127i 59.2 

0.05 3.965-0.241i 63.33 

0.25 13.175-1.021i 74.597 

0.45 26.99-2.19i 79.1046 

0.60 40-3.3 i 81.015 
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Fig.3 Scattering coefficients of linear pol and RX pol vs. incidence angles (vms=0.15). 
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Fig.4 Scattering coefficients of different vms vs. the incidence angles 
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Fig. 5 Surface roughness effects on the RR and LR scattering 
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Fig. 6 Correlation functions on the RR and LR scattering 


